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There are many optical experiments which measure linewidths, lifetimes, T\ times and T2 times, etc., but 
it is not always clear how these relate to each other in terms of the basic quantum mechanics. This paper 
gives a basic introduction to these experiments based on elementary time-dependent many-body perturbation 
theory in the random-phase approximation. Of primary interest is to be able to write down this theory for 
the general case of a nonequilibrium many-body system without assuming an equilibrium temperature. We 
do not want even to assume a perturbation around a quasiequilibrium system, but a system which is far from 
equilibrium, e.g. a shell in momentum space. This type of highly nonequilibrium system is often created in 
ultrafast laser experiments. 

It turns out that it is nontrivial to show a very basic experimental result, that optical lines are broadened 
in the shape of a Lorentzian which has width proportional to the instantaneous scattering rate, in the case 
of a nonequilibrium many-body system. The assumptions necessary to deduce this result in turn show the 
limits of its applicability. 



1 Time-Dependent Perturbation Theory 

The theory here is based on standard time-dependendent perturbation theory, e.g. as presented by Baym 
0. We review here the basic results given in Chapter 12 of Baym. 

We begin by assuming that the Hamiltonian of a quantum mechanical system is given by 

H = H + V, 

where Hq is the single-particle Hamiltonian and V is a scattering term which is small compared to Hq ■ The 
system of interest can be anything from free atoms or ions to carriers in a semiconductor or metal. We 
ignore all relativistic effects, however. 

We assume that the system has been prepared at time t = in the quantum mechanical state l^o)- At 
a later time t, the state is written as \ipt}- 

The Schrodinger equation gives the time evolution of the system as 

ih^ t ) = (H a + V)\^ t ) 

In the interaction representation, we define a new state \ip{t)} (t in parenthesis rather than subscript) given 

by 

m)) = e iHat/h \^t) 

and a new operator V(t) = e lH o t / h y e - lH ot/fi-^ j n this representation, the Schrodinger equation is rewritten 
as 

which has the advantage of not depending on Hq. 

Integrating this equation as a perturbation series, one can show that [jjj 

\m) = e- {t/R) i v(t)dt mo)) (i) 

= ( l + (l/t&) [ V(t')dt' + (l/ih) 2 [ dt' f dt"V(t')V(t") + ... ) |^(0)) 
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2 Connection to Time-Independent Perturbation Theory 

In general, we are interested in knowing the amount of depletion of the initial state after some time t. 
Abbreviating |^>(0)) = |0), we write 

(0\m) = (l + (l/ih){0\V\Q) I dt' 







{l/ihf f dt' f dt"J2(0\V\m)(m\V\0) 
Jo Jo m 



e (i/K)(E -E m )t' e (i/K)(E m -E )t" 



t t' t" 

+ (l/ihf [ dt' [ dt" I dt"'J2{0\V\n){n\V\m){m\V\0) 



JO JO 



Xe (i/K)(E -E n )t' e {i/H){E n -E m )t" e (i/H){E m -E )t'" + \ 



(2) 



Here we have inserted a sum over the complete set of states \m) (m\ = 1 which are assumed to be eigenstates 
of the Hamiltonian H . Then we have 

W(t)) = (i + (i/ih)t(o\v\o) + ^(i/ihft 2 \(o\v\o)\ 2 

+ (l/ih) T f dt' l(m|y|0)|2 e m)(E„-E m) t> Ui/H){E m -E )t' _ A 



~(l/ih)H 3 \{0\V\0)\ 3 



+2(i/ l nr(o\v\o)J2 fdt 1 f dt"-J- 

x ^ e (i/h)(E m -E )t' - + ...^ . 



|(mly|0)] 2 M/h)(E -E m )t' 

- E m + irj 



(3) 



In going from (|^) to (^), we have used some assumptions. The upper and lower bounds of each time integral 
give terms of the form 

f e (i/h)(E m -E )t' _ ^ _ 

The first, exponential term cancels the exponential term of the next time integral. The second term, -1, 
gives a vanishing term in the integration because when the last time integration over dt' is performed, it 
gives an integral of the form 

V(E m )dE m - l) 'f 1 ^ 1 ' 

V / — £j m ) 

Assuming T>(E m )\ (m\V |0) | 2 is continuous and slowly varying with E m , this integral vanishes because the 
leading term is odd, i.e. 

^ e i(E m -E )t/n i{E m - E )t/h 
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for [E m — Eq) ~ 0, and when (E m — Eq) ^> 0, the 1/(Eq — E m ) 2 term in the denominator and the fast 
oscillation of the e l ( E m~E )t/h ^ erm ^his integral. This is the "random-phase approximation" (RPA); 
essentially, it means that we ignore memory of the past and keep only terms from the upper bound of the 
time integrals. 

Note also that rigorously, to introduce the term ir] in the denominator, which allows us to treat the pole 
at E m = E , we must assume that V — Ve~ vt , where i] — > 0. Then we will have 



and 



JO 



e-'" e"" 4 dt'dt" = / e"" 4 -dt' 



(1 - 2r)t + 2r}t 2 + ... - 1) - 2(1 - rjt + \t 2 + ... - 1 



rj 2r] 2 rj 2 

If 



2r] 2 

etc., which gives us the same result as if we had simply done the integrals assuming V is time-independent 
and then inserted it] whereever it is needed to take care of a pole. 

The series in (0) has the same form as the series expansion of an exponential. By induction, we can write 



(0|V(t)) = exp 



.(w i (oivrn + S l(m|m|2 



E - E m + if] 



(4) 



This result (|j) is extremely useful. It is not an approximation, but is exact to all orders of t, within the 
limits of the RPA, for any time-independent V. The series inside the exponential is just the time-independent 
perturbation series for the energy correction due to the particle interactions (see, e.g. Baym |IJ, chapter 11), 
which is typically called "Rayleigh-Schrodinger" perturbation theory. 
By the Dirac formula, the second-order term is 

= A^-iT^ (5) 

where the second term is just the total scattering rate. A*4) = (0| V|0) is called the mean-field energy, A^ 2 ) 
is called the real self-energy and I^ 2 ) is called the imaginary self-energy. From (ji|) we have therefore 

(0\m) = e -(^)(A (1) +A< 2 >-*r< 2 > )t (6) 

The probability of being in the state |0), given by | (0\ip(t)) | 2 , decreases over time as e~( 2T / h ^ t due to out- 
scattering. Thus one can see the reason why the imaginary self-energy is associated with an out-scattering 
rate. 
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3 Absorption and Emission Line Shape in the Single-Particle Pic- 
ture 



We now imagine that we have a single particle prepared in a given state which can decay by two channels, 
namely by coupling to an external field (normally a photon field) and by scattering to other internal states. 
We write 

H = H + V + V 

where V is a smaller perturbation than V, which gives the interband photon-electron interaction. We will 
only be concerned about the first-order contributions of this term, i.e. we will not worry about renormal- 
ization of the electron states due to electron-photon interaction. We will allow renormalization due to V, 
however. Therefore we write 

mt)) = \iKo))+i f dt'V(t')m')), 

in Jq 

where IVK^')) i s assumed to include the evolution due to V given in ([!]). 

We are interested in the rate of emission of a photon with frequency u> and momentum k. We assume 
that the eigenstates of Hq are single-particle states, each with a given momentum. Therefore by momentum 
conservation, only one state with momentum k, which we will identify as |0), couples to the photon state 
|w,k). We write the rate of emission as the probability of being in the photon state |w,k) at t = oo, after 
starting in the initial state |0), divided by the total time spent in the initial state: 



1 



i( W ,k|^(oo))r 



T (u,k) j °°dt\{om))\ 2 

By the assumptions above, by momentum conservation we write (u>, k|V'|m) = (V')So_ m . Then we have 



|< W ,kh//(oo))| 2 = 



/ df y2(oj,k\e^ Hot y' e ^l^t \ m )(m\^{t')) 
Jo 



dt' (cj,k\V'\0)e luJt ' e - {l/K)Eat ' (0|V>(i')> 



= \(V)\ 2 

- KOI 2 
= KOI 2 



1 



Tiuo — E — iT 
1 



{huj - Ef + r 2 

where we have abbreviated E = Eq + + and T = I^ 2 ' from The normalization factor is 

2 



dt\(o\m)r = 



dt 



-(iAHAd'+A^-irt 2 ')* 



dte~ 2Tt / h 



h 

2T 



5 



This gives the total rate as 



m? ^ 2r ^ IT . 9 (7) 



r(w,k) lx " (fod - £) 2 + r 2 

which when T — > 0, is 

1 2"7T 

The interpretation of this result is that while the photon has definite momentum and energy, the particle 
emitting the photon has definite momentum but indefinite energy. The energy uncertainty comes from the 
Heisenberg uncertainty relation AEAt > h, that is, the shorter the time spent in state |0), the greater the 
energy uncertainty. 



4 Connection to the Boltzmann Equation in Many-Particle The- 
ory 

We have so far been vague about the nature of the interaction V. In the previous section we assumed that 
we had a set of single-particle states with a time-independent out-scattering term V. In general, however, 
the particle of interest will be one of many particles in a system, and V will be an interaction term between 
the particles. 

The theory of Sections [l] and || is completely general; that is, the state |0) can be taken as a many-particle 
state instead of a single-particle state. If it is taken as a many-particle state, however, then the result (^) 
is not as useful, since it involves the total self-energy of the system, while we are typically interested in the 
single-particle self-energy. 

We define the instantaneous many-particle state using creation and destruction operators a k and a k 
(Baym chapter 19) as 

\t 



(8) 



where n k gives the instantaneous occupation number of each state (this is called a "Fock" state). The 
"vacuum" state |vac) is the zero-particle state, which in the case of a solid means the ground state of 
the system. Note that we do not need to assume an equilibrium distribution of particles; we can use an 
instantaneous nonequilibrium distribution if we have that information. 

A typical interaction term is written in terms of the same creation and destruction operators, e.g. a 
two-body, number-conserving term, 

V = \ C/ k 1: k 2 ,k 3 ,k 4 a k4 a ] j < . 3 a k2 a kl (9) 

ki,k 2 ,k 3 

where the summation is not over k4 because it is implicitly assumed that momentum is conserved so that 
k4 = ki + k2 — k3. The interaction energy U is assumed to be symmetric on exchange of ki with k3 or k2 
with k4. 

We are concerned about the evolution of the single-particle state k, which is to say, the number of 
particles in state k as a function of time. The change in the number of particles in a time t is given by 

d(n k ) = (V*|nkhk> - <0|n k |0) 
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= (ij(t)\e lHot/h nue- lHot / h \m) ~ (0|»k|0) 

= (0| e (4/ri) / y(t)dt n ke - (i/?l) / y(t)dt |0) - (0|n k |0) 

= (0|e (4/ri) / y(t)dt [n k , e - (4/?l) / V(t)dt ]|0). (10) 

The operator ro k commutes with H$, by definition. If it commutes with V, then there is no change in (n k ) 
over time. We can resolve the commutator in ( |l0| ) by using the relations 

[n k , a k '] = -a k <5 kik ' 

[n k ,a k ,] = a k (5 kik/ (11) 

which are valid, surprisingly, for both boson and fermion creation and destruction operators. For a four- 
operator term in the interaction (ft), we have 



n ka k4 a ka a k2 a kl = al 3 a k2 a kl <5 k:k4 + a t4 a k3 a k2 a kl £) k:k3 
~ a k 4 a k 3 a k2 a ki<5k : k 2 - a k4 a k a k2 a kl (5 kjkl 



- a k 4 a k 3 ak 2 a ki n k 



Thus 

[n>k,V] = -^^2 (^032ia k a k3 ak 2 aki + f/302i a k3 a k a k2 a kl - E/32oia k3 a k2 a k a kl - %2ioa k3 a k2 akiak 



2 

ki,k 2 



^ ( U D ± ^b) (a^L^aka - al^l^au) (12) 



2 

ki,k 2 

where U o refers to the direct term and Ue to the exchange term, and the + sign is for bosons and the - sign 
is for fermions. (For hard-sphere scattering, i.e. s-wave scattering, U is a constant, which gives a factor of 4 
enhancement of the scattering cross section for bosons and is forbidden for fermions.) 
We first write out the series expansion, 

d(n k ) = (0| - (l/ih) J V(t')dt' + .. 

x [ (l/ih) [ dt'[n k ,V(t')] + (l/ih) 2 I dt' I dt"[n k ,V(t')V(t")} + ... I |0) 



Any terms in the right-hand scries multiplied by the leading "1" in the left-hand series vanish, since 
(0|[n k , A]\0) vanishes for any operator A. Note that V^(t) = V(t). The leading-order term is therefore 



ft \ / ct 

;2\ 



d(n k ) = (0\(l/h 2 )[ / dt'V(t')\[ / dt"[n k ,V(t")} |0) 



^(1/a 2 ) Qf dt > e m)(E -E m )A Qf dt/ , e (im(E m -E )t"^ ( \v\m)(m\[n k ,V]\0) 

/ M/h)(E Q -E m )t _ 1 \ / e -(i/h)(E -E m )t _ i\ 
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where the sum over states m is over all possible Fock states. The time-dependent factors are resolved using 
the identities 

t e imt _ i) U-ixt _ x \ sin 2 (xi/2) 
hm -± = lim 



S(x)2irt 



We therefore have 



d(n k ) = £<0|V|m)<m|[n k ,V]|0)^<J(^-.E TB ). (13) 

m 

Since the different Fock states are assumed orthonormal, the summation in V is eliminated because only 
the terms which couple |0) to |m) survive. Since a destruction operator ak acting to the right on a state 
with rik particles gives a factor y 7 ^, and a creation operator a k gives a factor \/l ± rt k , where the + sign is 
for bosons (stimulated emission) and the - sign is for fermions (Pauli exclusion), using ([l2]), ( |l3| ) therefore 
becomes 

= ^ \Ud ± ^g| 2 [^k 3 "-k 2 (l ± ^ki)(l ± "k) - » k n kl (l ± n k2 )(l ± n ka )] 

ki,k 2 

X(S(£; k + S kl -S k3 -E k3 ) (14) 

where we have assumed i is a very small quantity dt. This is the quantum Boltzmann equation, which gives 
the total rate of out-scattering for a state k. It is the same energy we would have written down if we had 
simply written the total scattering rate as the sum of the Fermi's Golden Rule rate, 

\( m \V\0)\ 2 5(E -E m ), 

n 

over all allowed processes of single-particle states |m). It is correct to order V 2 , but the limit t — > oo in ( |l3| ) 
implies that the time step dt must be long compared to the oscillation time %/{Eq — E m ), i.e. we are not 
concerned about behavior on time scales so short that the equivalent energy uncertainty is comparable to 
the typical collision energies. This is the random-phase approximation again. 

We can make one simplification of the Boltzmann equation by treating nk as a continuous variable, in 
which rtk is equal to its average value in an element of phase space d 3 k, and converting the summation to 
an integral, to get 

cJ(rlk) 27r 



/ d^ki d 3 k 2 \U D ± U E \ 2 [n k3 n k2 (l ± n kl )(l ± n k ) - n k n kl (l ± n k2 )(l ± n ka )] 
at n \(27r) J / J ' 

xS(E k + E Wl -E k3 -E k4 ). (15) 

We can make an additional simplification if we assume that the system is isotropic, i.e. the distribution 
function depends only on the magnitude of k and not on the direction. Then we can integrate analytically 
over all the angles to reduce this integral to just a double integral which can then be solved numerically for 
nonequilibrium isotropic distributions. This has been used to produce predictions for various nonequilibrium 
systems [|,|,|,|. 

These different rates clarify the difference between T\, T 2 etc. in experiments. The rate (14) gives the rate 



of depletion of a single quantum state, sometimes called the "dephasing" rate in optics. The time constant 
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for this decay is called "T2" in NMR terminology. The rate (|l5|) gives the rate of depletion of states with 
the same macroscopic characteristics as state k, and the time constant is called "Ti" in NMR terminology. 

We now have calculated the many-body equivalent of the formula (^J) for the decay rate. We can now 
use this in the calculation of line broadening. 

5 Line broadening in Many-Particle theory 

As in Section we are interested in the rate of emission of a photon with momentum k and frequency u> from 
a given initial state, via a weak particle-photon interaction in the presence of a particle-particle interaction. 
In the many-particle theory we will define the creation and destruction operators c k and Ck for photons; 
the total many-particle state is a product of Fock states of both photons and the interacting particles. As 
before, we write H = H + V + V, where V' is the particle-photon interaction term, which we assume has 
the form 

V = ]TP klCki a kl . 

ki 

Although this form is strictly unphysical, it is a good model for radiative transitions. Electrons are, of 
course, conserved, but if we restrict our attention to only one band, then radiative transitions which cause 
band-to-band transitions involve the destruction of an electron in one band. On the other hand, excitons are 
not conserved, so this notation allows us to treat electrons and excitons on equal footing as in the previous 
section. Also, we drop the hermitian conjugate term because we neglect the possibility of photon absorption; 
we assume there are no photons at t = 0. 

The number of photons emitted in a time t is then 

(iv k ) = (m\cW\m) 

= (0\e {l/h) I {v(t)+v ' i(t))dt [clc k ,e~ {l/h) J( y (*)+ v 'W) dt ]|0) (16) 

As in the previous section, we assume |0) is an exact Fock state which is not an eigenstate of the Hamiltonian. 
This will complicate the analysis (preventing us from using standard many-body formalism) but allows us 
to treat the case of a highly nonequilibrium system. As before, we write out the series, 

(JVk) = (o\(i-(i/%h)J {v{t') + v\t'))dt' 

+{l/ih) 2 jf* dt' J* dt" (V(t") + V'\t")) (V(t') + v'^t')) ...^ 
x ({1/iK) j\t'[clc^{V(t') + V\t'))] 

+ (1M) 2 jf* dt' jf* dt"[clc k , (V(f) + V'(t')) (V(t") + V\t")) } + ..^j |0) 

Since V is not number-conserving, only terms with products of the pair V'*V' survive. We restrict ourselves 
to only terms which are first order in V'^V. The photon number c^Ck commutes with V(t), while the 
relation ( |TT| ) implies that [c k Ck,V] = PkC^a^, which in turn implies that only terms with V'^ = ^ki k Ck 
survive. 
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At this point we can make an additional assumption in order to greatly simplify the calculation. This 
is to assume the "dilute" limit, which is that the occupation number of state k in |0) is not greater than 
one, and in-scattering to any particular state k is negligible compared to out-scattering. One can see that 
this will be the case because for any out-scattering term a^a^a^a^ for which the four momenta satisfy 
momentum conservation exactly, the probability of having both of states ki and k2 occupied simultaneously 
is very low. Note that in equilibrium, the average net scattering into and out of any state over a long time 
equals zero, but in this case the many-body state |0) has been prepared with n k = 1 at t = 0, and therefore 
the out-scattering at that time will greatly exceed the in-scattering. 

In terms like VV there are combinations like aj^a^ ak 2 a kiak- No double occupancy means that terms 
with ak<2k vanish. On the other hand, terms in V for which no k^ = k commute with ak and one can show 
that all such terms vanish. We will use Vk(t) to indicate the subset of terms in the summation of V(t) which 
do not commute with a^. We then get 



(JV k ) - (0||P k | 2 \-{l/ih) f aic^e-W^-^'dt' + {l/ihf f dt' f dt"V k (t")alc k 
V Jo Jo Jo 



e (i/h)(Hu>-E k )t' 



-{l/ihf f dt' f dt" f dt'''V k (t''')Vu(t'')alc k e-^ nu - E ^ t ' + 
Jo Jo Jo 

e t r t ft' 



x((l/ih)[ dt'cla^/^-^*' + (1/ih) 2 f dt' f dt'^la^^-^'v^t") 
\ Jo Jo Jo 

t t' t" \ 

+{l/ihf [ dt' [ dt" [ dt"'cla^ h ^- E ^ t 'v k {t")V ii (t'") + ... j |0) 
Jo Jo Jo J 

= |fk| 2 (o| (-(i/ih) e - {t,h){huJ - Ek)t ' dt 1 " 



x 1 1- (i/ih) f v k (t")dt" + (i/ih) 2 f f Vu(t")Vu{t"')dt"dt"' 
y Jo Jo Jo 



t' r t' r t" 

/f\ iii! i n /„.*\2 



x ( i + (i/ih) I v k (t")dt" + (i/ihf J J v k (t")v k (t"')dt"dt"' + ... I |o) 



h 2 



xn k 



e -(i/h)(hu-E k )t' dt ^ e (i/h) J*' v k (t")dt" 
(jf e {1/ih){hw ~ Ek)t ' dt'^ e~ {t/h) So W') dt " |o) (17) 



We now break this into two parts, inserting a complete set of states \ m )( m \ an d separating the term with 
|0)(0| from the rest. 



(AW - <-« 



(J e-WV^"-^' dt'(Q\e {i,h) SI Vk( - t '"> dt " 



10) 
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xn k 



' e (l/iH)(Hu-E k )t' dt n Q i e -(i/K)J* V k {t")dt" ^ 

f dt ' e -m)(.^-E k )t' f* df n ( |^ k (t")| 
/o Jo 

The second, cross term vanishes, having terms like 



m/0 



(18) 



, e iW2 *' - 1 



C1/2 



(Wl + W 2 )W2 UJ1UO2 



which vanish in the limit of the random-phase approximation, while the other factors are given by (|4|) - (|^). 
In the limit t — *■ 00, we then have 



(N k ) = ^>k) 
= Ia| 2 (n k > 



dt > e (i/h)(hu-E k )t' e -( J /r i )(A( 1 )+A( 2 >- l r< 2 >) 



(Rw - £) 2 + r 2 ' 

where we have used the abbreviated notation for E and T of Section and 



(19) 



= 7rJ2\(™\VK\0)\ 2 HEo-E m ) 



= *E 



\\ Y,( U D±U E )al 3 al 2 a Wl a w \0) 



k,,k 2 



<5(i?o — £ m ) 



= 7T J] |^D±C/B| 2 ^k 1 (l±nk 2 )(l±«k 3 )'5(Sk + £ ; k 1 -^ka --Efc.) 
ki,k 2 



(20) 



Note that the integral in ( p0| ) is the same as in (14), neglecting the in- flow terms, which are negligible in the 
dilute limit for a particular state k, even though the average in- flow to a region of phase space d 3 k around 
state k, found in (|15|), equals the out-flow from the same region in equilibrium. 
We can rewrite (|l9|) as 



dt 



= m 



, d(n k ) 



1 



dt (htu - e) 2 + r 2 
2 2r/a 



(nu - Ef + r 2 



(21) 



in which ([bfl ) is used to obtain d(rik)/dt = —2T/%. 

This all has been essentially a long justification of a very simple formula, which we already obtained 
in (0) using the single-particle picture. This calculation shows us the limits of validity of this formula, 
however. In particular, it assumes the dilute limit. In the gas of a high-density boson gas, the situation 
is much more complex, as shown by Shi, Verachaka and Griffin || for case of a quasiequilibrium boson 
gas. When state k is nondegenerate, however, this calculation shows that we can use the the interaction 
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term V = Yl \ Ud ± 3= ™2\/l ± "3 with generality in the single-particle picture for a particle in 

interacting with a bath that is in equilibrium or far from equilibrium. This type of interaction potential has 
been used to predict the Lorentzian line broadening of luminescence from various semiconductors ^| , known 
as "homogeneous" broadening (as opposed to "inhomogeneous" broadening, due to random fluctuations of 
the band gap.) 

An even simpler method, which is justified by this calculation, is to simply write H = Hq + V' where H' 
includes the self-energy correction to the single-particle energy due to scattering by V . In this case we would 
write e('/ fi ) H o*|0) = e^/ n ^ Bk+A k + A k - jr<2) )* anc [ then treat the interaction V in first-order perturbation 
theory. 

An important implication of the calculations here is that even in the case when there is substantial line 
broadening, the energy conserving delta- function in ( |l4|) should not be replaced by a Lorentzian. One can 
show that if one were to do this, total energy conservation of the system would be violated. 



6 The Connection to x in Four- Wave Mixing 

We have seen that the imaginary self energy gives the line broadening of optical transitions. Alternatively, 
as seen above, it also corresponds to the out-scattering rate. This rate can often be measured directly in 
time-resolved measurements. The two methods are complementary. If the rate is extremely fast, it may be 
too fast for time-resolved methods, but then it gives signficant line broadening. On the other hand, if the 
line broadening is too small for the experimental spectral resolution to pick up, it corresponds to slow decay 
rate in time-resolved measurements. 

The nonlinear susciptibility is defined in terms of Maxwell's wave equation in the presence of a polarization 
current </, 

8 2 E d 2 E dJ 

^=Moeo^+Mo^. (22) 

The polarization current depends on the electric field in a complex way, which is represented to third order 

as 

J = eo |(x£ + X (2) £ 2 +X (3) £ 3 ). (23) 

To calculate x' 3 ' from first principles, then, we need to calculate the current due to the oscillating dipole 
moment of the medium. If we recall that the J = eV = eP/m, we see that fundamentally we want to 
calculate the momentum (P) of the oscillator as a function of electric field. The value of can then be 
identified as the part of this which depends to the third order on the electric field E. 

For simplicity we will assume that the oscillator has two states which couple only weakly to other states. 
These two states nominally correspond to the valence and conduction electron bands in a semiconductor. 
We will label the ground state \v) and the excited state |c). We will also assume that the light field couples 
only the valence band to the excited band and has no effect on coupling within bands. 

Having justified its use in the previous section, we use the simple approach of replacing each state's energy 
with the appropriate renormalized energy including imaginary self-energy due to intraband interactions, 
subject to the same constraint of no multiple occupancy of states, and treat the optical transitions in first- 
order perturbation theory. Let the time-dependent electron state be |e t ). In the interaction representation 
we use |e(i)) = e lHat / h \et) . Then the oscillator momentum of interest is 

P(t) = (e t \P\e t ) = (e(t)\e lHat/h Pe- lHat/h \e(t)) (24) 
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= (e(i)|e^*/^K)(i|P^|i)(i|e-^*/ fi |e(i)) 

* 3 

= (e(t)\c}e l ^ +ir ^ t / h (c\P\v)(v\e(t)}e~ l{E -- tr ^ t/h 

+ {e{t)\v)e^ +lT ^ t ' h {v\P\c){c\e{t))e- l ^-^^ t ' n 
= 2^{v\P\c){e{t)\v){c\e{t))e^- E ^ t / h e- { - r - +T ^ t ' n , 

where (v\P\c) is the standard oscillator strength of the transition. 

The time dependence of the electronic state in response to the electric field is determined by the 
Schrodinger equation 

ih^-\e t ) = H\e t ) (25) 

H = H + V (26) 

V = -^—P-A (27) 
2mc y ' 

- - 4f ■ 

V' is the radiation term which connects the valence and conduction band states, and Ho is the total Hamil- 
tonian of the electron states. A and E are assumed to be classical fields. Let A(t) = A e~ lut . Then 

E{t) = (i/c)A Q uje- wt = E a e- wt (29) 

which implies 



In the interaction representation, 



v = e J—!k e -i"t (30) 



V'(t) = e iHot/h^O e -u,t e -iH t/h^ (31) 

v ; 2mw v ; 



Third-order time-dependent perturbation theory then gives 



|e(t)) = \e{t )) + ^j\t'V'{t')\e(t )) (32) 
dt' / dt"V'(t')V'(t")\e(t )) 

Ho J to 

ft' rt" 



1 



i- fdt' f dt" f dt"'V'(t')V'(t")V'(t"')\e(t Q )) 

lfl ) Jt Jto Jt 



(i 



where to - * — oo. 

Let |e(to)) = \v). Then since we assume that V does not couple intraband states, so that only (c|y|w) 
terms survive, we have 

(c\e(t)) = ±.j\t'(c\V>(t>)\v} (33) 
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(v\e(t)) 



1 + 



f dt' j dt" f dt"'{c\V'{t')\v)(v\V'{t")\c){c\V'{t"')\v) 

to J to Jto 

t pt' 

dt' / dt"(v\V'(t')\c)(c\V'(t")\v) 

to J to 



Therefore 



P{t) = 2^{v\P\c)e l ^- E ^ t ' n e- {T ^ +T ^ t ' h 
1 



- / dt'{c\V'{t')\v) 



(34) 



to 



l 



t ft' pt" 

dt' / dt" / dt" , (c\V'(t')\v)(v\V'(t")\c)(c\V'(t'")\v) 

to -J to J to 

' dt'{c\V'{t')\v)\ ( fdt> f dt"(c\V'(t')\v)(v\V'(t")\c) 

to J \J to J to , 



{my 

+0{V' 5 ) 

There is no second-order \^ by symmetry. We will concentrate on the first third-order term; the calculation 
of the second term gives similar results. 

Substituting in for V'{t), and performing the time integrals explicitly, we get 

Pit) = 2^\\(v\P\c)\ 2 e^- E ^/ h e-^ +r ^ h ^- f dt ' e i(E c -ir c )t'/h e -i u t> -i(E u +ir u )f /h 

+ \(v\P\ c )\^^- E ^ h e-^+ T ^ h ( i^lV-i- ['dt' f dt" f dt'" 

\2mujJ (ihf J tQ J to J to 



x e 



x e 



(E c -iT c )t'/h e -iut' -i{E v +iT v )t' /h\ f i{E v -iT v )t" /% -iuit" -i{E c +iT c )t" /h 



) 



i(E c -iT c )t"' /h -iuit'" -i(E v +iT v )t"' /h 



(35) 



23? 



\( v \P\ c )\* e i(E u -E e )t/h e -(r v +r e )t/ti 



eE 



i{E c -hui-Ev)t/h 



2muj (E c — E v — htu + i(T c + T v )) 
3 , „t ,t' 



(T v +T c )t/h 



+ \{ v \P\ c )^ e ^-^)t/n e -(T v+ T c) t,n ( 1 f dt , f df „ 
x ^ e i( E c-ir c )t' /h e -iu>t' e -i(E v +ir v )t' /h^ ^ e i(E v -ir v )t" /h e -iu>t" e -i(E c +iv c )t" /h^ 



i(E c -hu>-E v )t" /h 



{T v +T c )t"/h 



{E c -E v -hw + i{T c + T v )) 

eE e- iut 1 



(36) 



v\P\~^ r 



23? 
+\(v\P\c)\ 



2muj {E c - E v - hu + i{T c + T v )) 

iu)t \ 3 



4 / eE e 
\ 2mbj 



' {E c - E v - hu + i{T c + r„)) 



- I (J_ 

{E c -E v - 3hw + i3{T c + T v )) \2lj 

2(v c +v v )t/h 



(37) 
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The first term of the final result is the linear \ which goes into the index of refraction, and the second term 
is proportional to E 3 , i.e. it gives X ■ Here we have deduced only a frequency-tripling nonlinear effect 
with an ingoing and outgoing resonance, because we assumed only one input frequency. Of course, if we 
had written A(t) = Ai(e Ml * + e~ luJlt ) + A 2 (e 1 ^ 1 + e^ 2 *) + A 3 (e luJ3t + e~ luJ3t ) instead of A(t) = A er wt as 
we did, then we would get all the frequency mixing terms associated with third-order optics, e.g. four-wave 
mixing (FWM). 

As seen in (|35|), the prefactor e l { E v-E c )t/h e -{T v +r c )t/h m g^jj f fa e terms is exactly canceled out by a 
term with opposite sign in the exponent from the time integrals, so that the oscillation of the polarization 
is only at the driving frequency and the third-order harmonics. 

Suppose that instead of continuing on forever, the electric field is shut off at t = 0, i.e. 

««)-{V:>s <w 

Then for times t > 0, the integrals in ( |35|) are time-independent. The only time dependence in P(t) 
comes from the e l { E v-Ec)t/h e -{Tv+T c )t/h p re f ac tor. The polarization continues to oscillate at frequency 
(E v — E c ) /h, which is resonant at both u> and 3w (the sum frequency). If there is no scattering (damping) it 
will continue ringing forever. The imaginary part of the self energy gives the rate at which the phase-coherent 
polarization dies. This is why it can be called a "dephasing" rate. There is also the feature, seen in (|37|), 
that the polarization grows exponentially in time at twice the decay rate, which implies the counterintuitive 
prediction that a FWM signal will have a longer rise time if there is less dephasing. This has been seen 
experimentally ||, [l(J . 

In the case of delayed FWM, a third, "probe" wave follows two "pump" waves by some time delay. The 
strength of the fourth wave then gives a measure of the T 2 time. This is because the first two waves each 
create a polarization wave in the medium due to the linear term in (|35|), which decays according to the 
prefactor e -( r ^+ r <=) t / r \ These two waves then become contributions to E at later times, which mix with the 
probe wave via the the third-order term. Since the strengths of the pump polarization waves decay, so will 
the delayed FWM signal, with exactly the same time constant. 

We have seen then that the FWM dephasing rate and the line broadening are both controlled by the same 
imaginary self-energy, i.e. the T 2 time. This leads one to expect that the line width and the FWM dephasing 
rate should be inversely proportional. This is generally true, but can be incorrect when a continuum of 
states is excited, (e.g. by an ultrafast pulse with considerable Heisenberg energy uncertainty), which leads 
to interference between different oscillator frequencies, as seen in Ref. p"T|. 



7 Scaling Laws in Recent Experiments 

The electron-electron Coulomb scattering process represents a special case of the two-body scattering dis- 
cussed in Section ^, because the interaction cross section diverges in the case of zero momentum transfer, 
and this divergence must be removed self-consistently by a screening length which depends on the density 
and the instantaneous energy distribution of the electrons. A previous publication j^] presented the results 
of scaling laws for electron-electron scattering based on a Boltzmann integral calculation taking into ac- 
count the dependence of the self-consistent screening length on the density, using a method valid both for 
nonequilibrium and equilibrium electron plasmas. The main conclusions of that work are reviewed here. 

First, the proper integral for the electron-electron scattering depends on the type of experiment. There 
are three different integrals which relate to the experiments. The first is the "total" scattering rate, i.e. the 
imaginary self-energy calculated in (p(l) above. This integral determines the rate of decay of x^ 3 \ as shown 
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dephasing (T2) energy relaxation (Ti) 

low density: 

3d constant (l| n 1 / 2 BJlf, |l| 



2d constant n |17 



high density: 

3d n 1 / 3 jPj n 2 / 3 

2d n 1 / 2 [14] n E7 



Table 1: Density dependence of rates for electron-electron scattering predicted in Ref. Experimental 
confirmations are given in the references in the table. 

above. A second integral is given by ( |l5| ) weighted by the momentum exchanged by the electrons in a given 
scattering event. This typically controls the electron scattering rate determined in transport measurements. 



A third integral is given by (15) weighted by the energy exchanged by the electrons in a given scattering 
event. This determines the evolution of the carrier distribution function, which is typically recorded in time- 
resolved luminescence experiments. These three rates are approximately the same for hard-sphere scattering 
but are very different for Coulomb scattering. This means that the "relaxation time approximation," in 
which all scattering processes are assumed to be characterized by a single "relaxation time," is justified in 
the case of short-range interactions but breaks down completely for long-range Coulomb scattering. 

Second, the scaling laws as a function of temperature and density for these various integrals were deter- 
mined for both two dimensions and three dimensions. At that time, several experiments had been done which 
gave scaling laws consistent with the results of the calculations; recent experiments also fulfill the predictions 
of that theory. Table 1 gives a summary of the predicted density dependences and the experiments which 
have observed them. An important result is that this theory predicts that the dephasing rate is independent 
of density at low density, in basic agreement with recent experiments (Ref. [Q ; also seen indirectly in Ref. 
p8[.) The "high-density regime" is defined as the regime in which the classical screening length becomes 
comparable to the interparticle spacing. A full treatment of the scattering rate in this regime would require 
accounting for the quantum wavefunctions; however, the proper scaling law for the dephasing can be found 
by the simple assumption that the screening length is pinned at the interparticle spacing. This assumption 
also correctly gives the crossover from the high-density to low-density scaling regimes. 

A third main conclusion of Ref. || was that the screening length of the electrons scales with density 
in same way even for a highly noncquilibrium distribution. Quantum memory effects are important in 
determining the exact evolution of polarization, but the scaling laws should remain the same even when 
quantum memory effects are taken into account. 

Recent experiments showing a near-constant dephasing rate at low density have been interpreted in terms 
of optical phonon scattering |L2| . A full study of dephasing as a function of the excition photon energy and 
the temperature should distinguish between these two interpretations, since Ref. § also gave predictions 
for the temperature dependence of the scaling laws. 

8 Conclusions 

The iostropic quantum Boltzmann equation has been used to produce quantitative predictions for numerous 
experimental systems far from equilibrium B |l ||, || . These calculations show that the integrals which enter 
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into these calculations are the same as the integrals which are used for optical line broadening and four-wave 
mixing, for the cases of a Fermi or nondegenerate Bose gas. In the case of electron Coulomb scattering, great 
care must be taken to use the properly weighted integral for different experiments. 
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